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ABSTRACT 


With the increase in speed of moving vehicles there is a growing 
need for rapid, easy to interprete display of information concerning 
its dynamic state. 

With growing equipment stability a remedy can grow out of research 
and application of computer aided information presentation, especially 
the production of proper displays. While all information can be made 
to correspond to familiar schemes as close as desired new approaches to 
the general problem must Re sought to find an improved information 
presentation. This study is concerned with the application of computer 
generated graphics for one or more display consoles for aid in the 
operation of the vehicle. In the first chapter a thought experiment 
is performed and evaluated and a program structure for the operation of 
a flight vehicle is described. In the second chapter a major aspect of 
the general program is chosen for a demonstration: a 'flight vehicle 
to destination’ situation is subjected to an experimental investigation. 

A general purpose computer with a graphics terminal that has become 
available just recently has been chosen for the study. Thereby im- 
plicitly a measure of the Pasa lity to perform a real time on line job 


is obtained. 
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I. INTRODUCTION 


Flight vehicles including a human operator are made to respond to 
his controls. A standard way to present the information necessary for 
the operation of the particular flight vehicle has developed. Different 
vehicles (commerical aircraft, military aircraft, or even a submarine) 
have their favourable way to show the important data concerning the 
dynamic state of the vehicle. For the individual type there exist 
certain situations, in which control forces applied will be based on 
information which might probably be of relevance in this one situation 
only. 

The traditional way to present the necessary information is to 
enumerate all situations and provide a set of meters and indicators for 
each. Thus for n situations there will be n sets of meters. The design 
policy tries to use one particular meter in as many situations as possiblé 
and suitable. This procedure has several disadvantages; among aaa! 

(1) Supervision is made difficult by "clutter" effects the 
Operator cannot or only to a small degree overcome by selective techni- 
ques, if n is large. 

(2) The information presented will often require further imagi- 
native operator's work, i.e. the presentation is often far from Rates! 
due to indirectness. 

(3) The qualification of an operator must meet high requirements: 
combinatorial speed as well as spacial imaginative power. The last two 
are of special importance in aircraft operation. There will.always be 


long training periods. 


It has been claimed that a digital computer having a digital 
CRT as a peripheral device can take over almost all the combinatorial 
work in a complex situation, leaving the true supervisory task to the 
operator: 

(a) No ‘clutter’ effects need remain because the operator 
has a large range of display options; all non-relevant information is 
Suppressed. 

(b) The information presented can be made optimal by suitably 
transforming and evaluating measured data, avoiding any unwanted 
indirectness. 

(c) The requirements on spacial imaginative talent of the 
operator can be lowered by providing three-dimensional pictures of a 
simplified ('filtered') geometry, having no anticues. 

An application of display techniques in flight vehicle mission per- 
formance evaluation requires that the problem of presenting state infor- 
mation is restated and new ways to solve it are sought. Given the 
tremendous flexibility of a high speed digital computer it is reasonable 
to aim at more than a true copy of existing i Some information 
will always be presented in one dimensional form (e.g. a gyro compass). 
The capability of a digital CRT to allow information to be presented in 
meter form will be demonstrated in Chapter III. So far only limited use 
of the human capability to orient in three dimensions has been possible. 
A reason for this is the lack of convincing equipment. Contact analog 
devices ('Conalog'), see ref. 9, synthesize a three dimensional structure 
from measured data. Use of ‘Conalog’ ind te deers has been made with 
success. They are a remedy only in the absence of a true perspective 


picture of the environment. If true perspective pictures in any direction 


can be generated the ‘'Conalog' presentation will be possible as a 
by-product, although still very useful. 

Three dimensional displays combined with standard presentations are 
regarded but a first step in the direction of new information presenta- 
tion. (An example has been demonstrated in Chapter III.) Bute require- 
ment for a digitally controlled information presentation is the existence 
of a managing program. In Chapter II, fundamental aspects of such a 


program are investigated. 


IIT. PRELIMINARY CONSIDERATIONS 


A. MANAGEMENT PROGRAM STRUCTURE 
In the following the structure of a program (that could run on an 
arbitrary present day general purpose computer) is described in block 


diacnam. Lormam oce sive. 


B. DESCRIPTION OF BLOCK DIAGRAM 

Updating of all the variables that describe the state of the system 
and its references must be done at precise intervals. It also must be 
completed within time limits. 

A real time clock will initiate analog to digital conversion (ADC) 
at the precise times, following which the transforming mathematical 
operations will be carried out. Several levels of interrupts and pro- 
gram traps are necessary. levels are indicated by multiple arrows 
leading to the block called UPDATE. The main clock interval is deter- 
mined by the control requirements as well as the plant response. In- 
creased accuracy requires prediction which in turn requires higher 
sampling rates. Those problems are of no concern here as long as the 
computational speed of the computer is sufficient. 

For the storage and updating of the display matrix see under C. 
Under 'EVALUTE DISPLAY CONTROL’ the proper programming paths are set 
up that establish the connection between the display matrix (which 
contains all available data-variables as well as constants-in trans- 
formed form) and the storage routine, which transfers as the next step 
the selected data into the display buffer. Depending on the type of 


display chosen by the operator a projecting routine converting three 
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FIGURE 1 


dimensional coordinates into two dimensional coordinates of the focal 
plane will be inserted in the path from storage (display matrix) to 
the display buffer. The display clock determines the time when to 
refresh the vectors drawn ordsehamthe da#ajin theebufifiers 

‘SAMPLE ALL INPUT DEVICES’ performs two independent tasks: one type 
of inputs is meant as controls for the plant, while the other determines 
if and what change of display is desired. The latter part thereby pro- 
vides the parameters for ‘EVALUATE DISPLAY CONTROL’. After completion 
of sampling the computer returns to an idling state or whatever batch 
processing type of a job has been assigned (e.g. postflight evaluation 


preparation by storing certain data). 


C. DISPLAY POLICIES AND DATA STRUCTURES 
1. List Processing 

Most of the graphics terminals are capable of accepting data in 
forms of lists that are displayed item by item sequentially by means of 
interrupts. Certain bits in the data or certain register values con- 
trolling vector generation will determine a blanked (MOVE) or an un- 
blanked (DRAW) motion of the beam. Whenever a separate refreshing 
mechanism that is not using the vital registers of the central pro- 
cessing unit is available, list processing is the preferred way. (Note: 
The term ‘list processing’ as used herein should be distinguished from 
a common usage, i.e. referring to techniques of scanning tree organized 
symbol structures and using chained storage allocation. The lists 
referred to here are simply sequential strings of micro-commands.) The 
number of refreshing passes through the list must be appreciably greater 
than the number of changes made in the list. Hence it is preferable in 


terms of computing time minimization (or brightness improvement) to 
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generate a list which in turn is selene instead of building up a 
picture as it is made up of vectors, processing them sequentially as 
soon as they are computed. If no list can be established, it is likely 
that an improper information presentation has been chosen. (E.g. if a 
rotating vector is displayed and the angular velocity is so high that 
the whole circle area is illuminated, there is no longer any position 
information presented which can be eee monitored!) Proper display 
must conform with the plant dynamic control. In cases where noise is 
Superimposed on a vector, a smear out will result which will be resolved 
by the eye into an optical average. It is recommended, however, to 
insert an optimal estimator if there is enough computational time avail- 
able and analog averages are unsatisfactory. It violates the policy of 
presenting clutterfree graphics to allow os vectors to be displayed. 

From the above it follows that data for display are constants 
over an extended period of time (in terms of machine time). All avail- 
able data should already be stored in usable form before the real time 
process starts. Read only memory can be uSed if the data is constant 
for all time. A fraction of the data will vary permanently and requires 
central processor storage. These data are those used in the projecting 
routines which are in general functions of the flight vehicle's position 
in space. Permanently constant data are those comprising the static 
part of an indicator: scale, numbers, radial lines, titles, etc. Only 
one or at most a few numbers per meter need high speed read/write memory 
bacations. 

2. Compute and Draw Procedure 
A different situation exists if there is no way to refresh while 


update and sampling computation is carried out. In this case the display 
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clock will start the update operation at the end of a selected number of 
refreshing paths. The total managing work will be distributed over 
several clock pulses. Unless the number of vectors is small (as is the 
case in the programming project of Shaper IN list processing is still 
advantageous. A trade off must be made between the average computation 
time per vector and the time the program must be trapped while the 
vectors are drawn (times required to store and modify will add to the 


overall list processing time). This latter part will be further dis- 


cussed in Chapter III. 


D. PROPOSALS FOR DISPLAY IN A SPECIAL CASE 

Since the flight vehicle dynamics determine its most suitable 
information presentation for the purpose of operation, it is difficult 
to generalize when configurations for one or more display scopes are 
discussed. If the flight vehicle is identified as an aircraft and some 
departure from standard presentations is allowed, a typical situation 
can be imagined. See ref. 10. If the proposal for a display configu- 
ration leaves room for changes and amendments, it is justified for the 
purpose of clarity to look at a particular situation. Since flexibility 
is one of the attractions a digitally controlled CRT has, certain con- 
clusions which will hold for many easee: can be drawn. By describing 
what can be done in a particular case and by investigating the require- 
ments, a scheme is set up that can be used in a variety of cases. 

1. Scope Presentation 

Control of a flight vehicle can be based upon 
(1) no further information than a set of meters, 


(2) visual observation of the surroundings and a few 


indicators. 
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The first neglects the built-in capacity of the human being to 
orient in three dimensions if a proper picture is presented. The second 
has the disadvantage that the visual We ietens will not always exist 
and/or the time available is not sufficient to make out the given cues. 

It is proposed to construct a true perspective picture of a known 
environment with the help of a digital computer. The perspective view 
shall be the main information for the operator to base his control 
decision on. The knowledge of the enviroment could come from a topo- 
graphical map and could be aided by radar measurements and telemetry to 
ascertain relevance. If there is lack of tabular (mapped) information, 
the picture on the scope could be constructed from electronic measurements 
only. Combinations of 'Conalog’ elements and real perspective pictures 
should be left at the operator's choice, also whatever indicator is 
desired to observe closely. Since the number of scopes will always be 
limited some control over the projecting subroutine should be given to 
the operator: e.g. change from normal to wide angle projection (wide 
angle photography) to include a larger sector of the surroundings. 

2. Picture Generation 
The task of the computer would be to 

(1) generate a correlated set of three dimensional coordi- 
nates, i.e. evaluate a set of vectors in three dimensions, 

(2) project these three dimensions on a focal plane, 

(3) display the vectors obtained in focal plane coordinates. 
The picture generated by this method may not contain massive anticues 
and should suppress hidden lines if necessary. Since this picture is 
permanently updated when the flight vehicle moves, it provides an excel- 


lent method to aid the space orientation of an operator. See ref. 15. 


Jig) 


Programs to project a general set of vectors in three dimensions 
suppressing hidden lines exist (ref. 11). In Chapter III the implementa- 
tion of the general features (projection and updating) on the Adage 
graphics terminal is described, also a typical indicator is displayed 
simultaneously with the three dimensional Woe, projected in perspective. 

The second type of display under the control of the operator is 
the display of meters with one or a few vectors indicating the variable 
values. It seems to be unnecessary to simulate completely a standard 
indicator as long as the same information is presented. 

The configurations mentioned above can be combined easily. An 
interesting combination will result from Pseudo 3D included in real 
perspective projections. Proposals concerning this combination are made 


in Chaptenmeiti: 


PA? 
' 
4 1 


E. COMPUTATIONAL REQUIREMENTS , DATA REQUIREMENTS , LANGUAGE 

CONSIDERATIONS 

Data requirements for a three dimensional transformation and pro- 
jection need a study of its own and will not be discussed further. The 
storage requirements and operations on data comprising the two dimensional 
parameters for indicators are discussed in the following: 

See figure 2. Step (1) is a read operation which is done prior to 
the real time operation of the system. Read in will be a matrix which 
contains in fixed point form the parameters necessary to generate a 
standard type of indicator: 

- initial and final value of the variable 
- number of radial lines for round or horizontal and vertical 


lines for rectangular meters 
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- shape, round or rectangular 

- numbers to be shown 

- spacing of radial lines 

- linear or nonlinear spacing 

- title 

- configurations desired and offset for various cases, scale. 

A set of routines must be provided to work on these input values to 
generate a list of starting and end points of vectors for drawing opera- 
tions. In general every meter will have one set of vectors from which 
it is built. Only different offset values will be aeueG to the vectors 
to make them usable in various configurations. The offset can be 
applied digitally or by using hardware if available. Before offset 
values can be added, a multiplication by a scale factor is necessary. 
Scaling too can be done with hardware or software depending on the 
machine available. From the coordinate value x, stored permanently, the 
value x' for use in a chosen configuration is obtained by solving the 
equation x' = deltax + scale times x, where deltax is the offset value. 
Scale is the reduction factor of the normalized quantity x. 

The proper values for the parameter matrix must be found by trial 
and error after an intelligent initial guess. If a high level language, 
say FORTRAN, is used as the programming language it is difficult to 
regain all information if it is stored in one matrix only. Probably 
three or four matrices are necessary and a waste of storage will result. 
If, however, the prestart program can be written in assembly language, 
all information is most suitably packed into one matrix.. The words can 
easily be evaluated in the prestart phase, most information being in YES 


or NO form and hence requiring one bit only. 
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Under PRESTART all data are computed and stored in a static array. 
The managing program will reference to the columns of the static array. 
It will authorize storage of chosen static (constant) data into the 
display buffer, where also the dynamic part is stored. The variable 
part of the meter is obtained in UPDATE, where analog to digital 
operation is triggered and the value referenced to the static array 
column of the particular quantity. 

The composition of a display (called configuration) can be influenced 
from a keyboard. The operator will decide, if a prearranged configuration 
of several indicators together with perhaps a perspective graph or one 


peculiar one for his momentary needs is most suitable. 


F, MODIFICATIONS 

It is a matter of empirical design, which shape and visual impres- 
Sion for an indicator is most suitable. Human factor considerations 
are necessary. If no special copies of today's standard equipment are 
asked for (a special programming effort would be required), then it is 
possible to add and subtract easily different indicators for different 
variables by modifying the input parameter matrix. A general generation 
program will be possible for the static array if no special copies are 
asked for By the user. This might not be suitable for Tat ion of 
aircraft instrumentation. In this case the data of the static array 
have to be obtained in a different process, which could be the same 
digital CRT on which the meter is sketched with light pen and joystick. 
A general generation program, however, might be highly desired in an 


equipment testing situation where quick changes are often wanted. 


wy 


G. SUMMARY 
The center of vehicle mission performance evaluation with the aid of 
a digital computer and CRT is a managing program, around which subroutines 


are oriented: 


PRESTART, under which all static data are generated, 
- UPDATE, where all variables are updated, 


- SAMPLE INPUTS and EVALUATE; display control, 


transferring chosen data to display buffer and initiation of display. 
The problem of acquisition of information about the geographical environ- 
ment has not been tackled. In Chapter III such astructure is assumed to 


be given. 
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TfL. PROGRAMMING PROJECT 


A. COMPUTATION AND DRAWING OF A VECTOR 

This chapter describes a progrimming project which has been carried 
out to clarify and to provide with an experimental basis the considera- 
tions of Chapter II. The computer being used (ADAGE GRAPHICS TERMINAL) 


has a) word length of 30 bits. The basic calculation is carried out in 


a 
Fated 


a 15 bit format. Almost all special instructions are laid out for 
half-word operation. The basic machine cycle is 2 microseconds, See 
ref. 5. 

The capability of the AGT10 to have list display while other opera- 
tions are carried out is provided. Since the number of vectors has been 
kept low, a procedure has been chosen where every vector is displayed 
immediately after finishing its computation. No interrupts are used 
at meceane or at the end of 4 vector drawing. (See Appendix A on 
vector generation.) The program continues immediately after the vectors 
that control vector generation have been provided with the necessary 
values, control, scale, offset, intensity, xy-position. Refreshing and 
computation is done in one loop. This procedure allows change of data 
(due to rotation f.i.) per cycle, it is the most efficient way to 
generate a dynamic display. The minimum refreshing rate is 40 passes 
through the display list per second for no flicker; even if a picture 
is changing appreciably within a second still no ‘jumps’ will occur. 

If flickerless display can be established, it is possible to follow this 
procedure. For a dynamic display, however, it is not required to allow 


for changes at 40 times per second. A slower rate (less than 10 times 


WS 


per second) will do. If the vector displayed represents the outcome 
of an electronic measurement, a noisy vector will be drawn if the measured 
quantity contains low frequency noise (which is typical because only low 
frequency quantities can be displayed). 
No vector drawing, however, may begin before the old vector has 

been completed. This is nowserious restriction. Ifscomputat ion»is 
done between vector drawings, the vector generation can easily be con- 
trolled. To refresh at 40 times per second is a strong timing require- 
ment for the overall program: updating and sampling, etc., must be 
distributed over several cycles to avoid flicker. 

(1) All calculations must be carried out in minimum possible time. 
This excludes the use of floating point or double precision (30 bit) 
arithmetic. 

(2) Storage must be used to trade for gain in time. 

(3) No high level language can ne allowed for the managing pro- 
gram or the projection/rotation routines. It is not possible to use a 
high level language efficiently in conmection with word preparation for 
the various registers. Bit setting and checking is an often recurring 


Operation. A remedy could be the insertion of a small assembly language 


routine. 


B. GRAPH COMPLEXITY VERSUS VECTOR COMPUTATION-VECTOR DRAWING PROCEDURE 
The graph complexity is measured by the number of vectors to be 


drawn. Drawing a computed vector allows two approaches: 


1. Compute and Draw 


See flow diagram. 
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Take one xyz-set, project, MOVE; take next xyz-set, project, DRAW. 






Use mode bit control, no interrupts. (For further explanation see 
Appendix A.) Repeat this process as often as possible per second, at 
least 40 times. 

As long as the computational work requires more than 30 micro- 
seconds, the vector drawing time is independent of the vector length. 
A projection of a vector (two xyz-sets) with the help of CVTXY (see 
Appendix C, program listings) has been determined to take 220 micro- 
seconds. If no other computational work is done, then flicker will. 
occur as soon as a list requires more than 0.025 seconds. This is the 


@ase for 66 vectons. Thiaiseis thesnumber for the mast favowrable situ- 






project 





ation with no housekeeping of any kind required. If updating, etc., must 


‘be done the number of vectors will be reduced by the time required for 


4 
‘ 
\ 


housekeeping. 


2. List Display 


See flow diagram. 


Zh 






project 


Compute the list of xy positions in the focal plane, store them in a 


list, display the list stored as long as there is no change required. 
If a change is required, recompute the list. 

Roughly 200 microseconds are required for the housekeeping, 
bit-setting and storage into list per vector. Let the program be 
trapped for 30 microseconds (longest allowed vector), i.e. 60 ms for 
one vector. With 60 ms per vector, 417 vectors can be drawn. This 
assumes a list of vectors stored in proper form. 

3. Comparison 

For a static picture repeated computation allows a maximum of 
about 70 vectors while list display can process about 420 vectors. 
With list display it has been assumed that the program is trapped while 
vectors are drawn. If updating and rotations are necessary the number 
of vectors must be reduced. The rotational work is proportional to 
the total number of vectors. If .1 of the time must be reserved 63 
vectors are possible with approach (1), while 378 vectors could be 
drawn otherwise. To avoid pulsation, however, list computation must 
be carried out over several periods (only 60 vectors can be computed 
and stored per period) before a new list is ready for display. It is 


not possible to display a partially updated list because of distortion. 
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Hence, the reduction in the number of vectors drawable is approximately 
twice as large with procedure III.B.2. Still 336 vectors can be drawn. 

A remedy for caseIII.B.2 is possible by not requiring program 
trapping while vectors are drawn. In this case all updating and 
rotational computation can be done while drawing of vectors is accom- 
plished on a cycle stealing basis. All. 420 vectors can be drawn. 
Within slightly more than 7 periods a list is recomputable. Next the 


~“ 


results of sectionIII.B are tabulated. See Table I. 





Program Reduction 
Interrupt Trapped Max. No. in No. of 
EOV, EOL During of Vectors Vectors 
Vector Possible Required 
Drawing Li Wot . 
Comp. 
Repeated Me | No 
Computation 


ist 





Optional Yes 
Processing 
e N 
List ; Yes No 420 = ; 
Processing Reduction 
a = number of vectors in list 
£f = fraction of computation time required 


TABLE I 


C. PROJECTION OF A CHOSEN GEOMETRY 

The restrictions described above are silent directions for every 
programming effort that aims at implementation of a flight situation. 

A simple geometric structure in three dimensions has been constructed 
‘and serves as the data the managing program will assign to various sub- 


routines to operate on. See fig. 3. 
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The vector end points are sequentially transferred to a projecting 
routine that transforms three dimensional coordinates into planar 
coordinates of the focal plane. The routine that has been written is 
general in that it will project 3D coordinates onto the focal plane 
from any focal point selected. The distance of the focal point from 
the focal plane is allowed to vary. Also the focal vector may vary. 
Parameters required are: 

C (cL oe co): the coordinates of the focal point; 

alpha, angle of line of sight (focal vector) with the positive 
X-axis 5 

beta, angle of focal vector with y-axis; 

gamma, angle of focal vector with z-axis; 

d, distance between focal plane and observer. 


The projection program (called CVTXY) implements the following equations 


CaCl tome OSC) [d = DSTNC] 
x x 
=c +d cos 
qd, , B 
aE = < + d cos y 
KK = d| [(x-c,) cos Qa + ne) cos 8 + (c-c,) cos y] 
E= C. + KK (x-c,) 
= Come why Cy =C 
n ‘ (y e 
C = Cc. + KK (z-c,_) 


XX = [(§-4,) cos B - (n-a,) cos a}|sin y 
ebm [-(§-q,) cos y + (C-q,) cos a]|sin @ 
(for sin y = 0) 

YY = (C-q,)|sin 7 
(2) = (n-a,) [sin 8 (for sin y = 0) 


Equations (1) and (2) have not been implemented. 


ZS 


See figure 3. First the coordinates of Q are found. Q is the point of 
intersection of the focal plane with its normal through C. Computed then 
are the coordinates of P‘' which is obtained by intersecting the ray from 
P to C with the foacl plane. §€, n and € are then used to find the two- 
dimensional coordinates XX, YY in the focal plane. 

(0, 0) in the focal plane is given by Q. For further information 
see ref. 13. 

The projection program can be used by any user provided the following 
restrictions are met: 

given that the focal vector points to the origin, the linear 

extension of the geometric structure may not lie beyond 25% of the total 
word length used (15 bits, maximum number = 37777, implies a maximum 
linear extension of 10000.) . If the extensions become bigger, over- 
flow cannot be avoided. Overflow prevention operations have only 
partially been implemented to avoid time losses by a complicated com- 
putational structure. This structure will eventually lead to floating 


point arithmetic which had been excluded at the beginning. 


D. ROTATION 
The rotation program (called Dl) implements rotations (orthogonal 
transformations) around the y-axis and around the z-axis according to 


the following equations 
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( for rotation around the y-axis 
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For further information on orthogonal transormations, see ref. 7. 
Rotation around a general set of coordinate axes, defined implicitly 
by specifying Cl. Gees Co. alpha, beta, gamma, is demonstrated in a pro- 


gram called ROT. 


E. APPLICATION TO FLIGHT VEHICLE MISSION PERFORMANCE EVALUATION 

Projection and rotation are the basis of any perspective picture 
usage in performance evaluation. Since the programs are general from 
a mathematical point of view, they can be used to display any function 
of two variables. The part of the function to be projected must be 
broken up into a series of short vectors, depending in length the 
frequency content of the function. To have a satisfactory display of a 
three dimensional structure, it is in general necessary to solve the 
hidden line problem. This problem has not been tackled in this project 
for the following reasons: 

(a) The time requirement of a general hidden line suppressing 

routine will be beyond what can be tolerated in an application where 


refreshing must be done by the CPU. 
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(b) A simple geometrical structure of limited extension viewed 
from a large distance will not be destroyed by anticues, if the vari- 
ations in the z-direction are few and isolated. It follows that the 
routines provided in this project are not general enough to be used in 
connection with, say, engineering drawings which shall be viewed from 
all sides. If there is no requirement for a dynamic display, then a 
hidden line suppressing routine can be added to the existing routines. 
For the purpose of demonstrating the capability of a digitally con- 
trolled CRT in flight vehicle performance evaluation, they are not 
required. The following part of a flight simulation has been implemented 
(the pilot is imagined to be at the focal point): 

let the structure of the geometry be given as shown in fig. 3. 

Let it correspond to an airfield (one can think of a real 
display with a PSEUIO 3D). Upon approaching this field the eye of the 
pilot (he is allowed unshaded view into all directions) will sense the 
changed perspective with respect to the airfield as a result of his 
motion. Proceeding towards the airfield on the same flight level will 
change the perspective to a top view. Side translation will result 
in a rotation around the z-axis. The presentation of the field is 
permanent onicn means that the observer's eye is fixed to the destination 
and not following the course of the vehicle. It is a matter of a side 
translation to move the structure out of the scope, if it is desired to 
have the operator's eye restricted to a small viewing angle (solid angle) 
in the direction of the course of the vehicle. A program has been pro- 
vided that limits the x and y values when they are about to overflow 
(see Appendix). It has been demonstrated, that a permanent perspective 


display of the bearing of the destination has its attractions. 
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The sequence of perspective pictures that will result from a moving 
vehicle can be obtained in a computer application in two different ways: 
(1) One end of the focal vector is fixed on the object, the 
position in space to which the other end is connected is changing, hence 
the focal plane will follow a rotation and translation. Cc. oe co); 
alpha, beta, and gamma must be updated permanently. A straightforward 
solution will result. The computational work necessary to perform up- 
dating is lengthy and difficult to carry out in fixed point arithmetic. 
It requires evaluation of three dot products and evaluation of inverse 
cosines. 
(2) Instead of letting the aircraft move around the airfield, 
let the relative motion of the airfield around the aircraft be used as 
a basis for computing the perspective change of the airfield. Instead 
of allowing the airfield to move and being repositioned by coordinate 
transformation (translation), which is necessary in fixed point arith- 
metic, leave the airfield permanently as the origin of the focal vector. 
Use the computed change of perspective to rotate the geometrical 
structure. 
This simplification is only justified under certain conditions: 
(1) extensions of the geometric structure in view are small and 
confined to a certain region, 
(2) the (initial) S must be large compared to the highest objects 
in the structure. 
There are difficulties in providing magnification as a measure of 
changing distance encountered the solution of which is so far unsatifactory. 
Under the assumption that a satisfactory solution to the last problem 


can be found, the method has attractions: 


WS, 


(1) considerably less computational work, simple implementation, 
only one focal point and one value for each of C (cL a co)» alpha, 
beta, gamma is required; part of the projection program is repetitive 
and can be bypassed, 

(2) no coordinate translations are necessary. 

Beyond the advantages stated, the speed of the operation can later 
be DRO ed by attaching hardware which can accomplish rotations. 

In order to achieve the visual impression of an approach to an 
airfield, it is necessary to control the rotation of the structure to 
be projected according to the computed perspective changes. The simplest 
way to control rotations is to project along the absolute x-axis and 
rotate the geometry by the proper value initially for true geographical 
reference. From then on altering the vehicle's course and speed will 
have consistent effects on the perspective. 

The angles by which rotation must be carried out are obtained by a 
simple simulation routine. According to the motion in the xy-plane, 
the updated position is used to compute the angular change with respect 
to a reference line. See fig. 4. 

For a relative motion of the object from A to B the angle by which 
the object appears to be rotated is given by 3 > 

ies m/2 - arctan (yp/x,). 

When the managing routine calls the update routine, the value for 
the rotation is computed and stored. Before the vectors are displayed, 
an investigation is made if there is a need of an angular or horizontal 


rotation. If so, rotation is®carmied wuteeetore projection. 
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G. CRITICAL COMMENTS 
The overall managing program, given in block diagram form (fig. 5) 
can be regarded successful in the following aspects: 
(1) implementation of projection rotation and translation in 
a running process, 
(2) vector generation from the digital processor (DPR2) instead 
of using the display operator provided by the Adage Company, 
(3) effective demonstration of a section of a flight situation 
(see photographic sequence), 
(4) an initial step in the direction of limiting has been carried 
out (see Appendix C). 
The overall managing program leaves several aspects open for further 
Studies: 
(1) the flight situation demonstrated is still in an experimental 
state - 

(a) timing of managing is not optimized, 

(b) the minimum levels of rotation that provide a smoothly 
varying display have not been determined, 

(c) the aspects of final approach to the ‘air field' have 
not been investigated, a routine that magnifies the object proportion- 
ally to the distance and reduces the geometric extensions by software 
limiting need to be introduced, 

(2) a study is suggested for the generation of a set of vectors 
describing in a simple but proper way a chosen ‘target! area, which must 


be allowed to contain any three dimensional structure whatsoever. 
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H. PROPOSALS FOR APPLICATION OF THE TECHNIQUE 

Without further investigations an assessment concerning the general 
applicability of the technique in an arbitrary flight situation cannot 
be made. This is especially true with the problem of finding a set of 
vectors describing the geometry being unsolved. But there are certain 
cases in which the destination geometry is very simple; in a ship to 
ship guided missile problem there is no difficulty in defining the tar- 
get geometry. A lack of profile might be a problem here. A combination 
of PSEUDO 3D elements and real target description is no problem from a 
programming point of view and could be employed to ease the guiding 
problem greatly. When looking at the photographic results there is no 
way to tell which vectors belong to the real world and which are super- 
imposed. Once a geometric reference is set up a misinterpretation is 
impossible when geometric measurements are updated. 

There is one particular feature that should draw interest: by 
showing the perspective view of an object as seen from a moving obser- 
vation point (the view is no function of roll, pitch and yaw.) all the 
time, it is possible to train personnel in minimum time to bring a 
vehicle into a desired target area. This feature is effectively dem- 
onstrated (see photographs). The sequence of photos mentioned depicts 
a situation as shown in fig. 4., where the technique is described, that 
leads to the sequence. The photos show also the true course of the 
vehicle and its radar picture, superimposed in the same display. With- 
out familiarization with the significance of the different types of 
information combined, the display will be confusing. This suggests 
that certain types of familiar presentations might not be combined at 


all. They are shown here to demonstrate the dynamic capability of the 
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display, and to allow control of the perspective rotations in the frame 
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IV. CONCLUSIONS 


A programming project has been carried out to 
- study aspects of flight vehicle performance evaluation with 
the help of a digital CRT, 
- investigate within the realm of the application the capa- 
bility of the AGT10 to perform the assigned task. 

A simple technique for operation (or guiding) a flight vehicle has 
been described and implemented in an experimental form. Some possible 
fields of application have been pointed out, also proposals for further 
studies have been made. The capability of the machine to perform a 
certain amount of computational work in updating and sampling controls 
has been demonstrated. With increasing complexity of the geometric 
structure list processing has been found the way that leaves a maximum 
of computation time available for other work. Updating and refreshing 
in one loop will be below the rate required as soon as the number of 
vectors to be displayed becomes larger. Suggestions have been made to 
introduce hardware in the rotational work, this will allow for more 
computation time in a tight situation. 

Fixed point arithmetic must be employed if the refreshing rate 
necessary for flickerless display shall be maintained. With fixed 
point arithmetic there will always be a finite probability of over- 


flow if numerical restrictions are not met. 
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APPENDIX A 


VECTOR GENERATION FROM THE DIGITAL PROCESSOR 
There exists a program package called DSPLY, the so-called display 
Operator, provided by the Adage Company. This system of programs fur- 
nishes means of drawing vectors, points and simple structures the 
defining coordinates of which must be described in machine language 
format. The aid consists of providing statements which are resolved 
into sequences of machine language instructions, which guarantee proper 
timing in transferring values to the various registers involved in a 
vector drawing operation. By avoiding these problems some advantages 
are lost and new problems are created (large storage requirement because 
of its generality). Much understanding for the otherwise often meaning- 
less formats of the display operator is gained by studying the basic 
steps of vector generation directly from the DPR2. A dynamic display, 
however, requires the most efficient vector generation and requires 
more storage than the display operator leaves for use. 
The following registers (destinations) must be provided with a 

value to control every phase and aspect of vector drawing: 

Destination 10 (D10), upper 14 bits 

destination 11 (Dll), upper 15 bits 

destination 7 (D7), both half-words 

destination 6 (D6), lower 15 bits 

destination 5 (D5), both half-words. 
The sequence shown above is mandatory. If a value is already set, its 


repetition can be omitted, except the transfer to destination 5. The 
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most important registers are destination 10 and destination 5. The 
operating modes of the vector generator are controlled by the values in 
the upper 14 bits of destination 10. D5 is the destination register for 
the XY pair to describe one end of a vector. One can distinguish two 
major operating modes of the vector generation: flag bit control and 


mode bit control. 
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A. FLAG BIT CONTROL 

M4 = 1 directs the vector generator (AVG) to investigate bits 14 
and 29, if a vector shall be drawn (bit 29 = 1) or if a vector shall 
be moved (blanked motion; bit 29 = 0). A vector drawing requires a 
MOVE operation to the starting point and a DRAW operation to the vector 
end point. If the starting point is the end point of another vector, 
MOVE can be omitted. Bit 14 = 1 signals end of list (EOL), in which 
case the AVG is switched off after the last vector drawing. Various 
options exist under flag bit control mode: the details can be obtained 
from the "Subsystem Programming Specifications" for the AVG. In the 
following is described, how a list is most suitably processed in a 
majority of applications: 

If M3 = 1, at the end of a vector drawn an ‘end of vector' (EOV) 

interrupt occurs and control is transferred automatically to pivot 
location) //i/56en¥ Thess instruction 


MDIR'X V7.1 56 
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is executed. The instruction to be taken to the instruction register 
(IR) is the one stored in 77756. If previously 
MDO5 DATA 

(transfer from memory location DATA to destination 5) has been stored in 
77756, then the next instruction is transferred to destination 5 the 
effective address being now DATA + 1 - due to the indexing operation 
which incremented the contents of 77756. DATA is the starting address 
of the list of values that represent the beginning and the end of all 
vectors to be drawn. DATA might typically be a set of vectors to con- 
struct a subpicture. The transfer to D5 is stopped if a 1 bit in position 
14 is discovered signalling the end of the list. No option exists to 
rule out the end of list interrupt under flag bit control. 

When the interrupt occurs control is transferred to pivot location 
77757 and the instruction 

PSRew tal 7 D1 

is carried out. Whatever address is stored in 7/7757, the next instruction 
is taken from this address plus 1 ag accordance with the correct inter- 
pretation of JPSR. The interrupt Flip Flop remains set until a JUMP'I 
instruction is encountered which clears the interrupt. If the address 
stored in 77757 is another transfer sequence D100, DI1, ..., D5, the 
basic refreshing loop is set up. With the help of an end of list 
interrupt it is possible to sequence several independent lists of 
variable length. All lists together might form a picture of a complex 
SEEUCEUTE. 

Transfer of the list values is made on a cycle stealing basis. If 
all list values must be modified before displaying can start (e.g. 
rotation of a body), the basic refreshing loop as described above can 


easily be altered. 
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The following instructions set the proper bits to accomplish a 


vector generation from a list with EOV and EOL interrupts (also all 


other register values are included) and trigger processing: 


MD16'0 ADDRI1 /operating mode of AVG 
MD11 ADDR2 /scale 
MDO7 ADDR3 /offset 
MD#6 ADDR4 /intensity 
MDAR ADDR5 /store next list addr in 
Ly Jou 
ARMD tee /store next list addr in 
oi) for 
MDAR ADDR6 /store MDO5 DATA in 77756 
ARMD 11750 /store MDO5 DATA in 77756 
MDO5 DATA Jtrans fex first XYopair 
continue 
In locations ADDR1, ADDR2, ..., ADDR6 the following constants are 
stored: 
ADDR. : 66466 ! 4H 
ADDR2 ey 7 ae 
ADDR3 : 6 
ADDR4 : 2f66 
ADDR5 : 9 : HH TABLE 
ADDR6 : MDQ5 ' Hf DATA 
DATA xe meet 
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B. MODE BIT CONTROL 

Mode bit control is exerted whenever M4 = 0. Bits (14) and (29) are 
of no effect. A variety of options exists. One very useful for indivi- 
dual MOVE or DRAW operations is outlined here: 

No interrupts occur at the end of the vector drawing operation. 
Hence the program must assume what the vector lengths and the resultant 
drawing times will be. The interval between any two transfers to 
destination 10, 6, or 5 should be about 8 microseconds (the equivalent 
of 2 normal instructions). A transfer too early acts like an inhibit 
Operation and no vector at all is drawn. Provided the other registers 
have been assigned the desired value previously, the instruction 
MD19'O LOC 1 
will initiate a MOVE operation, where LOC 1 10400 : H 
and the instruction 
MD19'O LOC 2 
will initiate a DRAW operation, where LOC 2 : 14406 : Be 
Obviously a word must have been transferred to destination 5 between 
the above two instructions and a transfer is expected after the latter 
one. 

With the help of the AVG subsystem programming specifications it is 
easy to resolve the word at location LOC2 into its information content: 
Bits 2, 3, and 6 are set, i.e. M2 = Ml = 1, M4 = M3 = 0. M4 = 0 implies 
mode bit control, Ml = 1 implies draw a vector in two dimensions, bit 
6 = 1 implies scope I on. To utilize the various options provided, it 
is easy to make up the proper word. When changing from a MOVE to a 
DRAW and back to a MOVE it is necessary to interpose an AND operation 


to make sure that the proper bit (Ml in this case) is on or off 
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respectively. Otherwise all MOVE instructions are interpreted as DRAW 
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C. SAMPLING OF FUNCTION SWITCHES AND JOYSTICK 
1. Function Switches 
Very often joystick and function switches are used in connection 
with vector drawing operations. The function switches individually have 
a value of 1 as long as they are depressed only. If IC (9) (interface 
control register) is on, the upper 15 bits of source 5 are connected to 


function switches 1 to 15. An instruction sequence 


MDIC‘O'L 

46 ! 4H 

S5MD FSW 
FSW : d 


will initiate the transfer of the instantaneous state of S5 to the 
memory location FSW. With suitable AND operations the various bits can 
then be investigated. If IC (9) = 0, then the states of function switch 
16 and the foot pedals are connected to 55. Obviously the contents of 
source 5 can be transferred to any legal destination like AR, BR, 

2, wJoysizack 

A different procedure must be followed if it is desired to obtain 

the instantaneous value of the joystick potentiometers. There is no 
analog to digital converter available with the AGTIO. Instead a digital 
comparison value is converted to an analog one and a comparison is made. 
Source 4 bit 24 is on or off, depending on whether or not the differ- 
ence between the comparison value and the unknown quantity is greater 


or less than zero. The comparison value is sent to destination /, bits 
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Oto 14. A multiplexing operation (initiated by switching on bits 22, 
23, or 24 of destination 11) selects the desired variable for comparison. 
It is possible to compare the sum of two unknown voltages. A sequence 
of yes or no questions (10) leads to a digital value of 10 bits. When 
the proper voltage value has been found, it is normalized with the 
maximum number found at the joystick. A program called 'HEBEL' can be 
called. The sampled value is returned in the lower 15 bits of the 
A-register. When entering the subroutine, the A-register must contain 
the values used in multiplexing (20 for the x-voltage, 40 for the 
y-voltage). When the end values of the potentiometers are not set 

to their proper value, normalization will give a meaningless result. 

A check of proper adjustment can be made by assembling a program 

called ‘ATODE', pushing at execution joystick into any extreme position. 
A listing of the values in locations ATODE + 16, and ATODE + 17, should 


show the corresponding values. Programs HEBEL and ATODE are listed 


below. 
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APPENDIX B 


PROGRAMMING HINTS FOR ADAGE ASSEMBLY LANGUAGE 

Careful study of the programming instruction manual provided by the 
Adage Company will not provide adequate preparation for several 
minor difficulties. The following remarks are intended to save experi- 
mentation time when first working with the terminal. 

(1) Multiplication involves the extended arithmetic unit (EAU1). 
Tt uses more time than a normal instruction. Hence, if the instruction 
that succeeds uses the A-registor, too, a NOOP must immediately follow 
Cer multiply IAs eLluccilony 

The same holds for the divide and shift instructions. The latter 

ones require a NOOP, however, only in the case of the shift being more 
than a half-word. 

(2) MPYU LOC multiplies the lower 15 bits of the accumulator 
by the contents of the upper half of the word in location LOC; MPYL 
LOC multiplies by the lower half of the word in location LOC, assumes 
bits 1 to 14 equal to zero, but taking part in the multiplication 
Operation. 

(3) Local storage assignment when programming in assembly language: 
can be done with the help of a macro -- 


REPEAT 


ENDR 


The number replacing the three dots must be octal. 
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(4) Operating recommendations for AGT 10: 

Program execution. To execute a program on the AGT 10 it is 
necessary to load the program from magnetic tape first. A teletype 
statement to the monitor START ("program name", 0)! will direct the 
monitor to load the program from tape with the help of the loader. The 
tape unit is indicated by ""'. All external programs needed are loaded 
and linked, too. A statement 120 ! or GOTO (120) ! will execute the 
program loaded. 120 can be replaced by any number consituting a mean- 
ingful entry point. For further details see programming specifications 
for AMRMX (name of system). A load and go procedure is suitable for 
some programs which is initiated by 'program name'! 

Generation of a program using assembly language. At the end 
of an assembly, the object program is appended as the next file on the 
user's tape. The assembler will operate on the source language state- 
ments in the scratch pad (file 10, first 24 feet of tape immediately 
following the load point) of the user's tape. The following steps are 
necessary in the production of an object program: 

(a) write source program on the scratch pad, 

(b) load assembler from mag tape andexecute assembler (either 
START ("ADEPT", 0) and ADEPT!, or simply ADEPT! as a load and go). 

To ( ): To write text onto the scratch pad is done with the help 

of a program called EDIT. When EDIT is loaded and executed the state- 
ments written on the scratch pad can be controlled from teletype. A 
variety of instructions exist (see Manual on EDIT). After text ‘edition’ 
has been completed, a control statement 'B, return' will try to return 
control to the monitor (system). Since EDIT overlays part of the system, 


however, it is necessary to reload the system from tape. This is done 


ay 


from the control panel (OPC) with depression of the following sequence 


of buttons: 

HALT - RESET - RUN - clear buffer register and set 1001, = PUISE MS 
Not earlier than after the depression of PULSE 1 the loader (which is 
not overlayed by EDIT under legal circumstances) will begin a tape 
search for file 11, which contains the system. When file 11 is found, 
a teletype statement MO/DA/YR = is issued; it must be completed properly. 

Step (b) can now follow. 

To (b): At the end of the assembly, it is again necessary to reload 
the system from tape, since the storage requirement of ADEPT and 
further necessary programs are such that the system is overlayed. If 
by some unwanted or unexplainable action the loader has been destroyed, 
it must be reloaded from paper tape (Bootstrap). An indication of a 
destroyed loader is often tape motion before PULSE 1 has been depressed 


when trying to load from tape. The procedure then is: depress buttons 


RESET and BTSTR in that sequence, then start paper tape reader. 
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APPENDIX C 


PROGRAM LISTINGS AND SHORT DOCUMENTATION 

A. SUBROUTINE SCHB 

Implements a translation of a vector in one direction (y-direction). 
(If the x-direction is direction of translation, rotation by 90° is 
necessary.) When a limit is met, scissoring is performed by computing 
a new vector end point. 

The value by which translation shall be made must be stored in the 
accumulator before entering SCHB. 

SCH expects the first address of the four vector coordinates 


(XX1, YY1, XX2, YY2) in the buffer register (Pointer). 


ee et re re Sa te 


add | 
offset | 







mod, 
necess. 


offset; 


F 
pointer = 
- 


! 
err 


4 J 
ey ae a ia 


modify by | 
computing | 
one new set 
of coord. 





return 





See photograph on following page. Routine VSRS is used for the purpose 


Omedaemonstration. 
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B. SUBROUTINE ROT 

Computes rotational change of coordinates for rotations by ie 
around the x, y and z axis. Six entry points can be referenced to 
externally. The mathematical basis is shown in Chapter III. Input 
parameters are (common storage) PX, PY, PZ. New values are returned 


in the same locations. 


C. SUBROUTINE HEBEL 


Subroutine HEBEL has been described in Appendix A (joystick). 


D. SUBROUTINE CVTXY AND Dl 


Subroutines CVTXY and Dl have been described in Chapter III. 
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EXPUNIGE 

Title sor 

ENTRY ROTsDRLFT»DRGHT»DRXPS»sDRXNGsDRYPS»sDRYNG 

COMMON CXsCY¥sCZsPXsaPYsPZa9XXs YY sDSTNCs ALPHAs BET As GAMMAs SBETA 
COMMON FilsF2sTHETA 


RO te NOOP 

D1s NOOP 
MDAR®TL 
GB707!'H 
ARMD DTHETA 
JPSR SCOSA 
ARMD CTHETA 
MDAR DTHETA 
JPSR S$SINA 
ARMD STHETA 
MDARIN DTHETA 
JPSR $COSA 
ARMD CFI 
MDAR'N DTHETA 
JPSR S$SINA 
ARMD SFI 
MDAR'L 
NOOP > 
ARMD DRLFT+1 
ARMD DRGHT +1 
JUMP "I D1 

Wes NOOP 
MDAR®L 
QOTATYH | 
ARMD DMEGA 
JIPSR $COSA 
ARMD CEPS 
MDAR DMEGA 
weSR SSINA 
ARMD SEPS 
MDAR'N DMEGA 
JPSR S$COSA 
ARMD CKAPPA 
MDAR'N DMEGA 
JPSR- SSINA 
ARMD SKAPPA 
MDAR SL 
NUOP 
ARMD DRXPS+1 
ARMD DRXNG+1 
JUMP ‘I D2 

D3 NOOP 
MDAR'L 
QB7TO7T'H 
ARMD DLBDA 
oe Sk S$COSA 
ARMD CMU 
MDAR DLBDA 
JPSR SSINA 
ARMD SMU 
MDAR’N DLBDA 


56 
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CEPS 


[COMPUTED XPRIME »sYrPRimMe 


a7 


DRXNGS 


DRY BS: 


NOOP 
ARMD 
MDAR 
MPYU 
NOOP 
ARMD'N 
MDAR 
MPYU 
NOQP 
MDAE 
ARMD 'H 
MDAR 
MPYU 
NOOP 
MDAE 
ARMD *H 
JUMP ‘I 


NOOP 
ra adel 4 
MDAR 
MPYU 
NOOP 
ARMD 
MDAR 
MPYU 
NOOP 
ARMD'N 
MDAR 
MPru 
NOOP 
MDAE 
ARMD 'H 
MDAR 
MPYU 
NOOP 
MDAE 
ARMD °H 
JUMP °I 
NOOP 
JPSR 
MDAR 
Mir ye 
NOOP 
ARMD 
MDAR 
MPYU 
NOOP 
ARMD 
MDAR‘'N 
MP YU 
NOOP 
MDAE 
ARMD 'H 
MDAR 
MPYU 


PR1 
PY 
Sor Ss 


PR4 
PZ 
SEPS 


PR1 
mg 
7 
CEPS 


PR4 
Pe, 
BIRAPS 


De 
ey 
CKAPPA 


PRI 
PY 
SKAPPA 


PR4 
ae: 
SKAPPA 


PR1 

ray 

eZ, 
CKAPPA 


PR4 
PZ 
DRXNG 


D3 
PX 
CMU 


PR1 
PX 
SMU 


PR4 
fu, 
SMU 


PR1 
PX 
Paz 
CMU 
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DRYNGS 


DTHETAS 
STHETAS 
CTHETAS 
SFI: 
Crs 
Pies 
PRA s 
DMEGAS 
CEPSs: 
Sime ¢ 
CKAPPAS 
SKAPPAS 
DLBDAS 
CMUs 
SMU: 
CNUs 
Sates 


NOOP 
MDAE 
ARMD 


id 


JUMP 'I 


NOOP 
JPSR 
MDAR 
mPYU 
NOOP 
ARMD 
MDAR 
MPYU 
NOOP 
ARMD 
MDAR 
MPYU 
NOUP 
MDAE 
ARMD 
MDAR 
MPYU 
NOOP 
MDAE 
ARMD 
JUMP 


NOOP 
NOOP 
NOOP 
NOOP 
NOOP 
NOOP 
NOOP 
NOOP 
NOOP 
NOOP 
NOOP 
NOOP 
NOOP 
NOOP 
Na P 
NOOP 
NOOP 


MPYU=154040!H 
TERMINATE 


°H 


"H 
" 


PR4 
PZ 
DRYPS 


D3 
PX 
CNU 


rns J 
PX 
SNU 


PR4 
ec 
SNU 


PR1 
PX 
22 
CNU 


PR4 
i A 
DRYNG 
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EXPUNGE 


TITLE CVTXY 
ENTRY CVTXY 
CXsCYsCZaPXsPYsPZaXXaxYYaDSTNCsALPHAs BETAs GAMMAs SBETA 


COMMON 
CIMMON 
GVIXY: 


Nl 


FlsF2sTHETA 


NOOP 
MDAR 
Cie > 
JUMP 
MDAR 
Q 

ARMD 
MDAR 
JPSR 
ARMD 
MDAR 
JPSR 
ARMD 
MDAR 
JPSR 
ARMD 
MDAR 
JPSR 
ARMD 
MDAR 
ues 
ARMD 
MDAR 
MPYU 
NOOP 
ARRS 
MDAS 
ARMD 
MDAR 
MPYU 
NOOP 
ARRS 
MDAS 
ARMD 
MDAR 
MPYU 
NJOP 
ARRS 
MDAS 
ARMD 


F 
Ni+e2 
BEGIN 


a 
ALPHA 
SCOSA 
CALPH 
BETA 
SCOSA 
CBETA 
GAMMA 
SBCOSA 
CGAMA 
BETA 
SBSINA 
SBETA 
GAMMA 
SSINA 
SGAMA 
DSTNC 
CALPH 


iss 
CX 
QXx 
DSTNC 
CBETA 


17 
CY 
QY 
DSTNC 
CGAMA 


17 
CZ 
QZ 


CALPHA STORED IN BITS 0-14 


C COS ALPHA IN BITS @ TO 14 


{C COS BETA IN BITS @ TO 14 


C COS GAMMA IN BITS @ TO 14 


C SIN BETA IN BITS @ TO 14 


CE SIN GAMMA IN BITS @ TO 14 


EIN BITS "1S 0geo 


CLAST 6 SHEPS COMPUTE OX IN@iSeee 


CQOY IN BES! Serres 


C IN 15 TO 29 
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BEGIN: 


NIIP 
MDAR 
MPYU 
NOOP 
ARRS 
NOOP 
ARMD 
MDAR 
MPYU 
NOOP 
ARRS 
NOOP 
MDAS 
ARMD 
MDAR 
MPYU 
NQOQP 
ARRS 
NOUOP 
ARMD 
MDAR 
Mel ag) 
NOOP 
ARRS 
NJOP 
MDAS 
ARMD 
MDAR 
MPYU 
NOOP 
ARRS 
NOOP 
ARMD 
MDAR 
MPYU 
NOOP 
ARRS 
NOUP 
MDAS 
ARMD 
MDAS 
MDAS 
ARMD 
MDAR 
ARRS 
DIVU 
NOOP 
ARMD 


vin 


CX 
CALPH 


28 


XMCX 
eas 
CALPH 


20 


XMCX 
XMCX 
CY 
CBETA 


2a 


AOD ¢ 
Py 
CBETA 


20 


YMCY 
YMCY 
CZ 
CGAMA 


2 


ZMCZ 
PZ 
CGAMA 


29 


ZMCZ 
ZMCZ 
YMCY 
XMCX 
NENNR 
DSTNC 
1 
NENNR 


KKK 


Cis TO 29 
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ONE? 


MDAR 
MPYU 
NOOP 
ARMD 
MDAR 
MPYU 
NOOP 
MDAE 
MDAE 
ARMD 
MDAR 
MPYU 
NOOP 
ARMD 
MDAR 
MPYU 
NOOP 
MDAE 
MDAE 
ARMD 
MDAR 
MPYU 
NIOP 
ARMD 
MDAR 
MPYU 
NOOP 
MDAE 
MDAE 
ARMD 
MDAR 
MPYU 
NOOP 
ARMD 
MDAR 
ARRS 
MPYU 
NOOP 
MDAE 


"H 


"H 


*N 
*H 


ARMD > 


MDAR 
MPYU 
NOOP 
ARMD 
MDAR 
ARRS 
MPYU 
NOOP 
MDAE 
ARMD 
MDAE 
DIVU 
NOOP 
ARMD 
MDAR 


MDAE *H'N 


CZKKK 
Gz 
ZETA 
QY 
CALPH 


XXZLe 
ETA 
17 
CALPH 


XXZL2 
XXZLe 
QX 

CBETA 


XXZL1 
XI 
17 
CBETA 


XXZL1 
XXZL1 
XXZL2 
SGAMA 


XX 
ZETA 
QZ 


C@ TO 28 
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DIVU 


NOOP 

AED 

J IMP 
CALPH:s N) IP 
CBETA: NOP 
CGAMA: NOOP 
SGAMA:NOOP 
QXs NOOP 
QY: N JOP 
ie: NIOQP 
XACXs NOUP 
TMCY : NJOP 
Ze 7: NQOP 
AKKs NOOP 
NENNRs: NOOP 
a NOOP 
ETA: NOQP 
Zit TAS NOOP 
XXZL1I:s NOOP 
XXZL2@: NUVP 
CXKKKs NOOP 
CYykKKKs NOOP 
CZKKKs NOOP 
woke: NOOP 


ARLS=17599!H 
MPYIT=15440!H 
ARRS=174430!H 
MPYU=154930!H 
MPYL=15420!H 
DIVU=164060!H 
TERMINATE 


UL 


SGAMA 


ey. 
CVTXY 
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.XPUNGE 


WEE Vsks 
ENTRY VSRS 


VSRS3: 


RTURN: 


RAN ¢ 


Beeb: 


NOOP 
MDAR'L 
6800) 
ARMD 
MDAR'L 
TAG9B 
ARMD 
MDAR'L 
200029 
ARMD 
ARMD 
MDAR 
MDBR 
JPSR 
JPLS 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
MDBR 
JPSR 
JPLS 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR*'H 
MDAR'A'L 
T7777!H 
ARMD 'H 
MDAR- 
MDAR'A'L 
BYH77777 
MDAE'H 
ARMD 
BITCR 

Q 
MD10°O°L 
10409!H 
MD11°L 
30800!H 
MDO7°L 
i 
MD@6°L 


PARX1 


PARY 1 


PARX2 — 


PARY2 
DEY 
PTR1 
SSCHB 
SAMPL 
RXCST 
PARY1 
RYSML 
PARX1 
RXMOD 
PARY2 
RYBIG 
PARXe 
DELX 
PTR1 
S$SCHB 
SAMPL 
RY SML 
VXX1 
RXCST 
VYY1 
RYBIG 
VXXe 
RXMOD 
VYYe 
VXX1 


VXX1 
VYY1 


VXX1 
POSXY 


64 


SAMPL: 


Ms 


M2: 


Mais 
DY POS: 


DXPOS:3 


DXNEG3: 


PARXI13 
PARYI1: 
PARX23 
PARY23 
RXCSTs 


Gay 
MD9OS5 
MDAR‘H 
MDARTA'TL 
T7777!H 
ARMD°H 
MDAR 
MDAR* ASL 
Ogi tT 
MDAE*H 
ARMD 
BITCR 

i) 
MD19°O°L 
14409!H 
MD@S 
MOLG OL 
ZNGria| 
SSAR*°H'F 
MDAR'A'TL 
49000 
Jia Bae 
SaAak HF 
MDAR'A'TL 
20009 
JPLS 
SSAR*H'F 
MDAR‘TA'TL 
10000 
ees, 
SSAR*H*F 
MDAR'TA'TL 
40900 ~ 
aL S 
MDAR 
MDAS'F 
ARMD 
JUMP 
MDAR 
MDAS*N'F 


~ ARMD 


JUMP 
MDAR 
MDAS*F 
ARMD 
JUMP 
MDAR 
MDAS'N'F 
ARMD 
JUMP 


QeQooao a 


POSXY 
VXXe 


VXXe 
Verve 


VXXe 


BUSA Y 


POSXY 


DYPOS 


DYNEG 


DXPOS 
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RY Scie 
RXMODs: 
RY BiG 
Wares au Le 
VXX18 
Vieie > 
VXXKe38 
Dewar s 
DEE: 
Phos NOOUP!HPARX 1 
Pine: NOOR HRY SML 
POSAY: 6) 
BITCR=5001@!H 
MD35=25809!H 
MDY6=26800!H 
MDA7=27T99G!H 
MD19=30000!H 
MD11=31960!H 
TERMINATE 


Qe VQoaaaoad ® 
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EXPUNGE 
SCHB 


TITLE 
ENTRY 
SCHB: 


Flils 


Shes 


SCHB 
NOOP 
ARMD 


ARMD*N 


BRMD 


MDAR'TI 


ARMD 


MDAR‘X'I 


ARMD 


MDAR*X "I 


ARMD > 


MDAR*X ‘I 


ARMD 
MDAR 
JZES 
JUMP 


ARAR'TH'F 


JPAN | 
JUMP 


MDAR*H 


JPAN 


MDAR‘'H 


JPAN 


MDAE SH'N 


JPAN 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
ARMD 
JUMP 


MDAR*H 


JPAN 
JUMP 


MDAR‘H'N 
MDAE °H 


JPAN 
JUMP 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
ARMD 


MDAR'*N 
ARBR‘F 
MDAR‘N 


ARMD > 
BRMD 
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Gigi: 


G22. 


Gels 


G31s 


MDAR‘SN 
ARBR'IF 
MDAR'N 
ARMD 
BRMD 
MDAR'H 
JPAN 
MDAR'L 
37777!H 
MDAE*H'N 
MDAE'H'N 
JPAN 
MDAE °H'N 
JPAN 
JUMP 
ARARIN'F 
MDAE*H'N 
JPAN 
JUMP 
MDAE'L 
37777!H 
MDAE'H'N 
JPAN 
MDAR‘H 
JPAN 
MDAR'L 
Sittin on 
4DAE'H'N 
MDAESH'N 
JPAN 
MDAE ‘HN 
JPAN 
JUMP 
ARAR'N'F 
DIVL 
NOOP 
ARMD 
MDAR‘'L 
JUMP > 
ARMD 
JUMP 


ARAR'N‘'F 


DIVL 
NOOP 
ARMD 
MDAR'L 
NOOP 
ARMD 
JUMP 
ARAR'N'F 
MDAE 'H'N 
JPAN | 
JUMP 
MDAE'L 
37777!H 
MDAE'H'N 


XXCST 


XXMOD 
XXCST 
XXMOD 
YY SMis 
Gl 


YYSML 
DDY 
STF Sal 
DDY 
STP51 
Ge 


DDY 
G11 
Ge 


DDY 
STPog 
YYBIG 
G3 


YYBIG 
DDY 
Gel 
DDY 
G2e 
STF 4 


DDY 
FRI 
STF 36 
KIMP 
SLPS 
DDY 
FRI 
KIMP 
STF3 
DDY 


G31 
STF4 


DDY, 
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( 


G32: 


Sis: 


KIMP s$ 


JPAN 
JUMP 


ARAR'N'F 


PPVL 
NOOP 
ARMD 
MDAR'L 
JUMP 
ARMD 


MDAR‘H'N 


ARRS» 
ARMD 


MDAR'H'N 


ARRS 
MDAE‘'N 
ARMD 
ANAR‘'N 
ARMD 


MDAR'H'N 


ARRS 
ARMD 


MDAR‘TH'N 


ARRS» 
MDAE'N 
ARMD 
NOOP 
MDAR 
MPYU 
NOOP 
ARMD 
ANAR‘'N 
MDAE ‘N 
JPAN > 
MDAR 
ARRS 
DIVU 
NOOP 
ARMD *H 
My L 
NOOP 
ARAE 
ARMD 
MDAR‘H 
MDAE ‘N 
MDAE °N 
MDAE 'N 
ARMD *H 
JUMP. 
MDAR ‘N 
DIVU 
NOOP 
ARMD °H 
Miny lo 
NOOP 
MDAE ‘H 
MD AE 
ARMD ‘H 
JUMP 
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STE S6. 


He: 


STF 37s 


Sle 4: 


TORT: 


PSL Y: 


MDAR 
MPYL 
NOOP 


ARAR'H'F 


MPYU 
NOOP 
ARMD 
ANAR YN 
MDAE'N 
JPAN 
MDAR 
ARRS 
DIVU 
NOOP 
ARMD 
MDAR 
MDAS'N 
MDAS'N 
ARMD 
JUMP 
MDAR'N 
DIVU 
NOOP 
MDAS 
ARMD 
MDAR'F 
ARMD 
MDAR 
ARAS'F 
MDAS 
ARMD 
JUMP 
NOP 
MDAR 
ARAS TF 
MDAS 
4ARMD 
MDAR 
ARAS'F 
MDAS 
ARMD 
MDAR*N 
ARBR'IF 


MDAR'N 


ARMD 
BRMD 
MDAR'N 
ARBR'IF 
MDAR'N 
ARMD 
BRMD 
JUMP 
MDAR‘H 
JPAN | 
MDAR‘'H 
JPAN 


XCH1 
XXMOD 
XCH1 
XCH1 
XXMOD 
STFS37 
NUMR 
DZAE 


XXMOD 
XXMOD 
Co fa Ha | 
YYBIG 
DDY 


YYSML 
YYSML 
TORT 


DDY 


YYBIG 
YYBIG 
DDY 


YYSML 
Yeuorit. 
YYSML 


YYBIG 
YYSML 
YYBIG 
XXCST 


XXMOD 
XXCST 
XXMOD 
DORT 
YY2 
Al 
YY1 
Alll 
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a 


AYQ1: 


Ev iolss 


0 le es 


Sime: 


Pt }/s 


B2e: 


MDAESHTN 
JPAN 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
ARMD 
JUMP 
MDAR‘H 
JPAN 
JUMP 
MDARTHIN 
MDAE °H 
JPAN 

SO) es 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR 
ARMD 
MDAR‘'H 
JPAN 
MDAR ‘L 
37777!H 
MDAE THN 
MDAE THN 
JPAN 
MDAE'H'N 
JPAN 
JUMP 
ARARIN'F 
MDAE SHIN 
JP AN 
JUMP 
MDAE TL 
1 MME ee 
MDAE*H'N 
JPAN 
MDAR *H 
JPAN 
MDAR'L 
37777!H 
MDAE°H'N 
MDAE*H'N 
JPAN 
MDAE *H'N 
JPAN 
JUMP 
ARAR'N’F 
DIVL 


Ye 
All} 
Ye 
YYBIG 
Ye 
YYone 
XX i 
XXMOD 
XXe 
AXCST 
StRe2 
YY} 
All 
AG) 
NCe 
YY 1 
All} 
AQ} 
Ye 

Y (eG 
YY 1 
Yeas. 
XX2 
XXMOD 
XX} 
XXCST 
Y “SME 
Bi 


YYSML 
DEL Y 

STP51 
DEE 

Sao) 
Be 


DELY 
Bll 
Be 
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NOQP 


ARAD FRI 
MDAR'L 
JUMP STP36 
ARMD CIMP 
JUMP STP3 
B21: ARAR'N'F 
DIVL DELY 
NOuUP 
ARMD FRI 
MDAR'L 
NOOP 
ARMD CIMP 
JUMP STP3 
83: ARAR'N'F 
MDAE HEN DELY 
JPAN B31 
JUMP STP4 
331: MDAE*L 
3 Tama | 
MDAE HN DELY 
JPAN B32 
JUMP STP4 
332: ARARIN'F 
DIVL DELY 
NJOP 
ARMD FR1 
MDAR'L 
JUMP STP36 
ARMD CIMP 
STP3: MDAR'H YY1 
ARRS 1 
ARMD DABVE 
MDAR*H YY2 
ARRS 1 
MDAE'N DABVE 
ARMD DZAE 
ANAR'N DZAE 
ARMD DNOM 
MDAR'H XX] 
ARRS 1 
ARMD DBLOW 
MDAR°H XX2 3 
ARRS 1 
MDAE'N DBLOW 
ARMD > DNR 
CIMP: NUOP 
MDAR DELY 
MPYU DNR 
NOOP 
ARMD NUMR 
ANAR'N NUMR 
MDAE 'N DNOM 
JPAN Ca 
MDAR NUMR 
ARRS 1 
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Go 


SIE EROS: f 


STP37: 


OTVU 
NQIUP 
ARMD 
MPYL 
NOOP 
ARAER 
ARMD 
MDAR 
MDAE 
MDAE 
MDAE 
ARMD 
JUMP 
MDAR 
DIVU 
No) 
ARMD 
MPYL 
NOVP 
MDAE 
“DAE 
ARMD 
ei = 
MDAK 
mie yi: 
New) i 
ARAR 
MPYU 
NOOP 
ARMD 
ANAR 
MDAE 
JPAN 
MDAR 
ARRS 
DIVE 
NQOP 
ARMD 
MDAR 
MDAS 
MDAS 
ARMD 
JUMP 
MDAR 
DIVU 
NIQOP 
MDAS 
ARMD 
MDAR 
ARMD 
MDAR 
ARAS 
MDAS 
ARMD 
JUMP 


al 


"N 
oN 


"F 


XCH2 
XXMIQD 
XCH1 
XCH1 
XCHe 
XXMOD 
Sheed 
NUMR 
DZAE 


XCH1 
FR 1 


AXMUD 
XCH1 
XXMOD 
5 PPM 
PRI 
DEE 


DNR 


NMR 
NUMR 
DNOM 
C2 
NUMR 
| 
DZAE 


XCH1 
XXMOD 
XCH1 
XCH1 
XAMOD 
Sieot 
NUMR 
DZAE 


XXMOD 
XXM0D 
Sarg ¢ 
YYBIG 
DELY 


YYSML 
asm L 
DORT 


13 


STP45 


DORT 8 


STP Sas 


JUTCKS 


XX1°8 
XX28 
XXMODs: 
KXXCSTs: 
YYTs 
Vie: 
YYBIG: 
YY Svile: 
DARL e = 
Hare 
DABVEs: 
DZAEs: 
DNOMS 
DBLOWs 
DNRS 
NUMRs 
NR1s 
DAT1: 
DAT2: 
DAT33 
DAT43 


NOOP 
MDAR 
ARAS'F 
MDAS 
ARMD 
MDAR 
ARAS'F 
MDAS 
ARMD 
MDAR 


ARMD'X'I 


MDAR 


ARMD'X'I 


MDAR 


ARMD'X'I 


MDAR 


ARMD'X'I 


MDAR*‘F 
simi * 1 
MDAR'F 
JUMP TT 
MDAR 


ARMD'X*I 


MDAR 


ARMD'X "I 


MDAR 


ARMD'X'I 


MDAR 


ARMD'X'I 


MDAR‘'F 
JUMP "I 


QW Wer Feaa05 geowv@teoeoeveaecdcsd 


DELY 


YYBIG 
YYBIG 
Dis ley, 


YYSML 
YYSML 
XXCST 
SAVE 
MY Sine 
SAVE 
XXMOD 
SAVE 
YYBr'G 
SAVE 
0 
SCHB 
1 
SCHB 
XX 1 
SAVE 
v Yad 
SAVE 
XX2 
SAVE 
rY2 
SAVE 


SCHB 


74 





XKCH13 y 
ACHe: 9 
Dupes J 
SAVE: J 
$APYIJS15499TH 
meyvL=15427084 
DIVJ=164)949!H 
DIVL=164e9!4 
x=] 144977 
TerR4INATE 
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=KXPUNGE 

TITLE ATODE 

ENTRY ATODE 

ATODE: NOOP 
MDAR'L 
2YVIVGDI!HODBG2A 
JPSR 
MDAR'TATL 
BOBOOIHT7777 
ARMD 
MDAR'L 
29AAO!'HAGDABAD 
Jit sie 
MDARTATL 
VAABIVHTT7T77 
ARMD 
JUMP 'I 

XVAL?: 0 

Ne Y 

TERMINATE 


SHE BEL 


XVAL 


SH a 


YVAL 
ATODE 


76 





Eee aN. Ge 


CDRS HEBEL 
Oa a) dial =) De De 


Heese: 


Gus 


G) 
ee 


DONS 


eciT $ 


TAFELS: 


NJOP 
ARMD 
MDAR TL 
NJVPPHTAFEL 
ARMD 
MDO7'H 
MD1ITL 
TITTIT!HHI9DII9 
MDARTL 

‘) 

ARLI'F 
MDAR'L 
OQ'H77766 
ARMD 
ADB 
BRAR IF 
SALA 
JUMP 
MDAS'X'I 
ARBRIF 
ARAT'F*H 
4ADAR'X 
Jicg les 
JUMP 
MDAS'N'I 
MDAS'TXK'I 
ARBR'F 
ARQT'F'H 
MDAR °X 
limes 
JUMP 
BRARTHIN'TF 
DIVI 
NOOP 
som 
O'HTAFEL 
Seltabk 
290990 
14094 
4393 
2004 
19903 

400 

249 

100 

43 


GI+1 


TAFEL 
MF SCL 


Grr 
MF SCL 


44 
DOWN 
TAFEL 


CTR 
G1 
OUT 
TAFEL 
TAFEL 


CER 
G1 
QUT 
13499 


aijowmy ob 


MFSCL:3: 43000 
DIVI=16449!H 
MD11=31000!H 
ARAT=27100!H 
On11=31109°H 
MD10=309990!H 
MDOA7=2TBBB!IH 
Cans NOOQP 
TERMINATE 


10. 


Li: 


ie 


ie 
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